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PREFACE 


The work described in this report was performed by the Earth and Space 
Sciences Division of the Jet Propulsion Laboratory. 

This report is Chapter 3 of Application and Fundamentals of Turbulence, 
to be published by Plenum Press of London and New York. The report was pm- 
sented as a Short Course Lecture at the University of Tennessee Space Institute, 
Tullahoma, Tennessee, on January II, 1977. 
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ABSTRACT 


A review is given ot the application of linear stability theory ti- 
the problem of bounda r v- layer transition in incompressible Tlw. The theorv 
is put into a form suitable for three-dimensional boundary layers; both the 
temporal and spatial theories are examined; and a generalized Gas ter rcla- 
tion for three-dimensional boundary layers is derived. Numerical examp ie'^ 
include the stability characteristics of Falkner-Skan boundary layers, the 
accuracy of the two-dimensional Caster relation tor these boundary layers, 
and the magnitude and direction of the group velocity for oblique waves in 
the Blasius boundary layer. A review is given of the available experiments 
which bear on the validity of stability theory and its relation to transit ion. 

The final section is devoted to the application of stability theory to tnn- 

n ^ 

sition prediction. Lienmann's method, the e method, and the modified -• 
method, where n is related to the external disturbance level, are all dis- 
cussed. A different type of method, called the amplitude method, is described 
in which the wide-hand disturbance amplitude in the boundarv layer is esti- 
mated from stability theory and an interaction relation for the initial ampli- 
tude density of the most unstable frequency. This me,t;l^a is applied to thr 
effect of freestream turbulence on the transition of Fi Ikner-Sknn boundary 
layers . 
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I. HISTORICAL BACKGROUND 

The earliest explanation for the appearance of turbulence was that the 

laminar flow becomes unstable, and the linear stability theory was first 

developed to explore this possibility. A series of papers by Rayleigh^^^ 

provided many notable results concerning the instability of invlscid flows, 

such as inflectional instability, but little progress was made toward the 

original goal. Viscosity was commonly thought to act only to stabilize the 

( 2 ) 

flow, but in 1921 Prandtl showed that viscosity can also be destabilizing. 

It was this discovery that finally provided a mechanism for the instability 

of boundary layers in zero and favorable pressure gradients which are stable 

to purely inviscld disturbances. However, it was not until some years later 
(3) 

that Tollmien worked out a complete theory of boundary-layer stability, 

and for the first time computed a meaningful critical Reynolds number (Re^^), 

i.e. the lowest Reynolds number at which instability appears. Any expectation 

that instability and transition to turbulence are synonymous in boundary 

layers was dashed by the low value of for the flat-plate boundary layer. 

Tollmlen's calculation gave a value of 420 for the critical displacement- 

* * * 4 

thickness Reynolds number, which is equivalent to ® UjX /v * 6 x 10 . 

Even in the high turbulence level wind tunnels of that time, transition 

was found between Re^ = 3.5 x 10^ and 1 x 10^. 

In what can be considered the earliest application of linear stability 

(4) 

theory to transition prediction, Schlichtlng calculated the amplitude 
ratio A/A^ of the most amplified frequency as a function of Reynolds number 
for a flat-plate boundary layer, and found that this quantity had values between 
five and nine at the observed transition Reynolds numbers. Outside of 
Germany, the stability theory received little acceptance because of failure to 
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observe the predicted waves, mathematical difficulties, and also the feeling 
that a linear theory could not have much to say about the origin of turbulence 
which Is Inherently non-llnear. The experiment of Schubauer and Skramstad^^^ 
completely revised this opinion, and unequivocally demonstrated the existence 
of Instability waves In a boundary layer, their connection with transition, and 
the quantitative description of their behavior by the theory of Tollmlen and 
Schllchtlng. This experiment made an enormous impact nt the time of Its 
publication, and by its very completeness seemed to answer most questions 
concerning the linear theory. To a large extent, subsequent experimental 


work on transition went in other directions, and the possibility that linear 
theory can be quantitatively related to transition has not received a decisive 
experimental test. On the other hand, it is generally accepted that f.low 
parameters such as pressure gradient, suction and heat transfer qualitatively 
affect transition in the manner predicted by stability theory, and in 
particular that a flow predicted to be stable by the theory should remain 
laminar. This expectation has often oeen deceived. A good introduction 
into the complexity of transition and the difficulties Involved in trying 
to arrive at a rational approach to its prediction can be found in a report 
by Morkovin.^^^ 

Investigators in Germany applied the sf^o^lity theory to boundary layers 
with pressure gradients and suction, and this work is summarized in Schlichting's 
book.^^^ We may make particular mention of Pretsch's^^^ work, as he provided 
the only large body of numerical results for exact boundary- layer solutions 
before the advent of the computer age by calculating the stability characteris- 
tics of the Falkner-Skan family of velocity profiles. The mather.atics of the 

(9) 

asymptotic theory were put on a firmer foundation by Lin, and this work has 
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been continued by Reid.^^^^ When in about 1960 the digital computer reached 

a stage of development permitting the direct numerical solution of the primary 

differential equations, the linear theory was extended to many more boundary- 

layer flows: three-dimensional boundary layers (Brown^^^^); f ree-convection 

(12) Q3) 

boundary layers (Kurtz and Crandall and Nachtsheim ); compressible 
boundary layers (Brown^^^^ and Mack^^^^); boundary layers on compliant walls 

n 6^ 

(Landahl and Kaplan' }; a recomputation of the Falkner-Skan flows (Wazzan, 

Okamura, and Smith^^^^); a quasi-steady calculation of unsteady boundary layers 

( 18) ^ 

(Obremski, Morkovin, and Landahl '); and heated-wall water boundary layers 

(19) 

(Wazzan, Okamura, and Smith )* 

It will be the main purpose of this chapt r to explain in detail the 

use of linear stability theory as a means of transition prediction. Enough 

of the theory is presented in Section II to make it clear how the essential 

quantity, the amplitude ratio is obtained. The use of temporal and 

spatial amplification theories is discussed, and a numerical procedure presented 

which allows eigenvalues to be calculated to arbitrarily high Reynolds numbers. 

A few numerical examples arr given, and in Section III the available experiments 

bearing on stability theory are examined. Finally, in Section IV the appllca- 

9 

tion of linear theory to traneition prediction is taken up, and the e , 

9 

modified e and amplitude methods arc discussed and applied to the effect of 
freestream turbulence on transition. It must be emphasized that the subject 
matter is restricted to incompressible boundary layers along an impermeable 
surface of zero curvature in the absence of body forces. Gbrtler instability, 
free shear flows, and stratified, rotating and compressible flows, are all 
excluded* 

Before proceeding further, it is well to mention some general references. 
These are review articles on stability theory bv Schlicht ing^^^\ Shen^^^\ 
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Stuart^^^^ and Reid^^^^; and books by and Betchov and Criminale. 

Schlichting*s^^^ book on boundary-layer theory contains two chapters on stability 

f 261 

and transition* and Monln and Yaglom’s'* '' book on turbulence contains a slnglf' 

(27) 

lengthy chapter on the same subject* as does the book by White on viscous 

(28) (29) 

flew theory. Reviews of transition have been given by Dryden* Tani. 

Morkovln^^^ and Reshotko. An extensive discussion of both stability theory 

and transition* not all at high speeds in spite of the title* may be found in 

(31) 

the Morkovin-Mack recorded lectures . 

II. STABILITY THEORY 
A. Formulation of Eigenvalue Problem 
1. D erivation of Equations 

The stability theory starts with the time dependent Navier-Stokes 
equations* not the boundary-layer equations. We will restrict ourselves to 
the flow of a single incompressible fluid on a surface of negligible curvature. 
This simplification eliminates many possible sources of instability, but 
preserves the twdl^rtilch are essential to an understanding of the subject: 
inflectional and viscous instability. The Navier-Stokes equations for a 
viscous incompressible fluid in Cartesian coordinates are 



3u 

dx 


U + v 


9u 

ay 


+ w 


au 

az 


i iE 

p ax 


+ vV u. 


( 1 ) 


p 3y 


2 

vV V, 


( 2 ) 


3w , 3w . 3w , 3w 


192 +, 

p 3z 


■jV^w , 


(3) 


3u . ^ ^ 

3x 3y 3z 


0 . 


(4) 
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The x-axls Is In the direction of the freestream velocity, the y axis is 
noirmal to the surface, and the z axis is normal to the x y axes. The 
velocities u, v, w are in the x, y, z directions, respectively. The density 
Is p, the pressure p, and the kinematic viscosity coefficient v = p/p. 'Hie 
first three equations are the x, z, z momentum equations; the fourth equation 
the equation of continuity. 

All flow quantities are divided Into a steady mean-flow term and an 
unsteady fluctuation term. A typical term Is 

q(x,y,z,t) = Q(x,y,z) + q'(x,y,z,t). (5) 


The mean-flow terms satisfy the boundary-layer equations. When expressions 
similar to (5) are substituted Into (l)-(4) for all flow variables. e mean- 
flow terms dropped which are negligible by the boundary-layer ec^ttons, the 
mean boundary-layer equations subtracted out, and rhe nonlinear terms neglected, 
a much simplified, but still too col, plicated, system of equati<^fl results. The 
additional assumption of locall> parall-' flow. 


U = U(y), W » W(y) , V = 0, 


( 6 ) 


reduces the equations sufficiently so that upon the introduction of a sinusoidal 
disturbance they become ordinary differential equations. 

The parallel-flow equations are, in dimensionless form. 


m' 



3u' 

^ „ 3u' , dU 

+ W - — + V -r~ ” 

.iiE' 

3t 

3x 

3z dy 

p 3x 

?v* 

3v* 

1 w ' 

_ i iz' 

3t ^ 

3x 


p 3y 


(7) 

( 8 ) 
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aw* 

at 
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3w* 

ax 


+ H 


aw* 

az 


+ V* 


dy 


1 3e’ 

p az 


X 1 -r2 ’ 

+ -7w 


> 


(9) 


au* ^ ^ 

ax ay 3z 


0. 


( 10 ) 


Fron now until Section IV dimensional quantities will be denoted by asterisks. 

* * / .- 

The velocity scale is the freestream velocity U, , the length scale is L (left 

^ 2 

unspecified for the present), and the pressure scale . The Reynolds 
nuokber is 


* * * , V 

R « L /v . (11) 

These equations, which are the basis of almost all stability investigations, 

are exact for the flow in a channel or for Couette flow, but are only an 

approximation for boundary-layer flows. If the multiple scale or two-timing 

(32) 

TOthod first applied to this problem by Bouthier and later by Saric and 
(33) 

Nayfeh is used, the above equations appear as the zeroth approximation. 

The next approximation takes the growth of the boundary layer into account. 
Although this method may be desirable for refined calculations, it is not 
needed for an elementary presentation of the subject and will not be pursued 
here. 

2. Introduction of Sinusoidal Disturbances 

The final form of the differential equations, where the coefficients 
are functions only of y, and x, z, t appear only as derivatives, suggests 
the following type of disturbance: 

(u* ,v* ,w' ,p') * (f,<Ji,h,Tr) exp [ i(ax + Bz - (ot)) . (12) 


-fr 
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Here, f(y), ^(y), h(y), ir(y) are the complex aiiq>licude functions of the 
disturbance flow variables u', v', w', p*; a and 0 are the dimensionless wave 

^ ifi ^ ^ 

numbers 2xL /X and 2xL /X , where X and X are the wavelengths in the x and 

A Z X z 

* * * 

X directions, respectively, and u> is the dimensionless frequency (o L /\J^. For 
the moment, a, 0, u may be either real or complex. 

When (12) is substituted into (7)-(10), the following equations for 
the amplitude function are obtained: 


i(aU + 0W-<i)) f + U'<> » - law + i I f" - (o^ + 0^) fj , (13) 

i(aU + 0 W-U 1 ) (> = - w* + i U" - (a^ + 0^)(^1 , (U) 

l(aU + 0W-to) h + W’4» «* - 10W + i (h” - (a^ + 0^) hj , (15) 

i(of + 0h) + 4.' - 0. (16) 


The primes now refer to differentiation with respect to y. The boundary 
conditions are that the no-slip condition applies at the wall, 

f(0) = 0, <K0) « 0, h(0) = 0, (17J 

and that the disturbances go to zero (or are at least bounded) as y -»■ <», 

f(y) 0, <Ky) -*■ 0, h(y) -*■ 0 as y -»■ ®. (18) 

Since all of the boundary conditions are homogeneous, it can be expected that 
solutions to (13)- (16) will exist only for particular combinations of R, a, B, 


- 7 - 
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and (i!. Consequently, we have an eigenvalue problem, and the primary task Is to 
evaluate the eigenvalue relation 

g(o, 6, u, R) « 0 (19) 

for one parameter In terms of the others. 

3. Orr-Somnerfeld Equation 

If (13) and (IS) are combined to form af -f Sh, this conblnatlon can be 
eliminated by (16), and, after differentiation, ir’ can be eliminated by (14) 
to give 


Iv 2 2 2 2 2 

♦ ^ - 2(a + 6 ) (>•• + (a + 6^)^ ^ 

= iR|(aU + BW-w) { (|i" - (a^ + B^) <)1 

- (aU" + BM*’) (20) 

When W ■ 0, this equation reduces to the equation obtained by Squire, and 
when In addition B " 0, to 

.Iv . 2^" ^ 4^ 

(j> - 2a () + o iji 

» IR [ (oU-u) (<t)’’-a^4>) - aU" <^>] . (21) 

This is the Orr-Sommerfeld equation and is the basis for most of the work done 
in incomprcr'slble stability theory. 

The Orr-Sommerfeld equation is a fourth-order equation and applies to a 
two-dimensional boundary layer. However, we can observe that (20), which is 
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the equation for a three-dimensional disturbance in a three-dimensional boundary 
layer, is also a fourth-order equation. This fact can be exploited to relate 
more general cases to (21). We will Illustrate this possibility with Squire's 
equation. With W » 0, the transformation 


2 2-2 . 
a + 6 aa ,Ra^Ra,(iKXBt()a 


( 22 ) 


reduces (20) to 


- 2a^(J»"+ - IR ( (oU-w) ((|>"-a^(|>) 

- , (23) 

which is identical to (21), but in the transformed variables. Since U is 
unchanged, it is evident that if a and 3 are real, a three-dimensional 
stability problem at Reynolds number R has been reduced to a two-dimensional 
problem at the lower Reynolds number R. This is the celebrated Squire 
theorem which states that in a two-dimensional boundary layer with real wave 
numbers, instability appears first for a two-dimensional dfsturbance. Further- 
more, if b) = a)(a, R) has been determined for a given U(y), then u = (o(a, 3, R) 
is immediately known from (22) . 

However, if a and 3 are complex or the boundary layer is three dimensional, 
the ut’llty of the above transformation is lost. In the first instance, R is 
complex, and in the second the boundary-layer profile is not invariant under 
the transformation. In both of these cases, there is little point in proceeding 
beyond (20) . The important conclusion is that in stability problems governed by 
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(7)-(10), the di termination of the eigenvalues only requires the solution of 
a fourth-order equation. 

4. System of First-Order Equations 

Since there are numerous stability problems that cannot be reduced to 
a fourth-orJer >ystem, a more flexible approach is to abandon the Orr-Sommerfeld 
equation altog'ither and work in terms of a system of first-order equations. This 
approach can h^: illustrated with (7)-(10) although it doesn't reveal its full 
advantage until the eigenvalue problem is of higher than fourth order. 

Let 


of + ^ “ af, ah - Bf ■ a!i . 


(24) 


By adding and subtracting (7) and (9), the following replacement equations can 
be formed for these two linear combinations: 


;(aU + 3W-oj) af + (aU' + BW) <|. 

= - iiT(a^ + B^) + ^ + 6^) if) , (25) 

i(oU + 6W-(o) ^h + (aW - BU') 

- i (ih - (a^ + B^) oh) . (26) 


With, 


Z, = of, Z„ * of, Z- ■ 4i, Z, = V, 


Z j » ah , Zg “ ah ' , 


(27) 


-10 
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the four equations (25) » (8) , (26) , (10) can be written 


Zj - Z2 . < 28 ) 

22 • + 1R(«U + IW-«)I 2j^ 

+ (oD' + (su > A Zj + 1 ( 0 ^ + B^) R 2^, (29) 

Z 3 - - 1 Z^. (30) 

2 2 

2 ^ » - I 2^ - ( IfoU + BH-w) + ° ^ 1 Zy (31) 

Z 5 • Zg. (32) 

Zg - (aW’-BU’) R Zj 

+ la^ + + IRUU + BH-w)l Z^ . (33) 


The fact tha** first four of taese equations do not contain or Z^ 

confinns that the eigenvalues can be obtained from a fourth*order system even 
though we are really dealing nere with a sixth-order system* It is only the 
determination of all of the eigenfunctions that requires the solution of the 
full sixth-order system. This formulation is applicable when a and & are 
complex and to three-dimensional boundary layers. The quantity a has been 
introduced only to connect with other formulations, but has not been assigned a 
meaning. When a and 6 are real, a is obviously the wave number in the direction 
of wave propagation. 


- 11 - 
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Ir the freestrean, (28)- (33) have constant coefficients and thus solutions 
cf the ion 


Z^^^(y) « exp (X^y) . 


(3A) 


The characteristic values occur in pairs, and are easily determined to be 




n y 1/2 

A ■ + B (»U, + BW,-u)l 

3,4 1 1 


S.6 “ ^3.4 * 


(35) 

(36) 

(37) 


where Uj^ and are the freestream values of U(y) and W(y). Only the upper 
sign satisfies the boundary conditions at y -*■ «. The characteristic functions 
for Xj^ are 


A<J) - - 1(«' + . 


,<1) 


2 2 
+ B^) , 


*<» - 1 . 


» i(aUj + pW^-uj)/(a^ + , 


- 0 , A^J^ - 0 . 


(38) 

(39) 

(40) 

(41) 

(42) 
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For real a, 6 and <d this solution Is the linearized potential flow over a wavy 

2 2 1/2 

wall moving in the oirection of the wave nuid)er vector with velocity u/(a -I- B ) 

It can be called the Inviscid solution, althoi^h this Interpretation is valid 
only in the freestream. The characteristic functions for are 



ila^ + + iR(«Uj^ + , 


(43) 



(o^ + + iR(oU^ + , 


(44) 



(45) 





(46) 


This solution represents a viscous wave and can be called the viscous solution. 
The third solution is another viscous solution, and is 


A 


(5) 

1 


- 0 , 





(47) 


- 1 . (48) 

« - (o^ + B^ + lR(aHj^ + BMj-w)l^^^ (49) 


These three linearly Independent solutions are the key to the numerical method 
of obtaining eigenvalues as they provide the initial conditions of the numerical 
Integration. 
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B« Temporal and Spatial Theory 

If a, B and o) are all real, the disturbance propagates through the 
parallel mean flow with constant rms amplitude. If a and B are real, and o) 
is complex, the amplitude will change with time; if a and B are complex, 
and 0 ) is real, the amplitude will change with x and z. The former case is 
referred to as the temporal amplification theory, the latter as the spatial 
amplification theory. If all three quantities are complex, the disturbance 
will grow in both space and time. The original, and for many years the only, 
form of the theory was the temporal theory. However, in a steady mean flow 
the amplitude at a fixed point is Independent of time and it changes only 
with distance. The spatial theory gives this amplitude change in a more 
direct manner than does the temporal theory. 

1. Temporal Amplification Theory 

With 0 ) s + lb), and ex and B real, the disturbance can be written 
r i 

q*(x,y,2,t) “ q(y) exp (w^t) exp [i(ax + Bz - u)^t)l . (50) 

The magnitude of the wave number vector is 

a = , (51) 

and the angle between the direction of a and the x axis is 

» tan ^ (8/or) (52) 

The phase velocity, or the velocity with which the crests move normal to 
themselves, is 


-14 
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c . “ M^/oi . (53) 

pn r 

If A represents the magnitude of q* some particular y, say the y for which 
|q* I Is a maximum, then it follows from (50) that 

(l/A)(dA/dt) = . (54) 

He can identify as the temporal amplification rate. Obviously A could have 
been chosen at any y, and (54) would be the same. It is this property that 
enables us to talk about the "amplitude" of an instability wave in the same 
manner as the amplitude of a water wave even though this amplitude is a function 
of y. He may distinguish three possible cases: 


«! < 0 

damped disturbances » 


01^ « 0 

neutral disturbances* 

155) 

> 0 

amplified disturbances* 



The complex frequency may be written 

u * ac * a(Cj, + i c^) . (.56) 

The real part of c is equal to the phase velocity c^^, and ac^ is the temporal 
amplification rate. The quantity c appears frequently in the literature of 
stability theory. However, it cannot be used in the spatial theory, and since 
wave theory usxially employs oi and lo, with the phase velocity being Introduced 
as needed, we will adopt the same procedure. 
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2. Spatial Amplification Theory 

In the spatial theory, (ii Is real and the wave numbers In the x and z 
directions are complex. With 

a - + 1 Oj, B - + 1 Sj , (57) 

we can write the disturbance In the form 


q*(x,y,z,t) “ q(y) exp [-(Oj^x + B^^z)! exp [ l(a^x + B^z - wt)] . (58) 


By analogy with the temporal theory, we may define a real wave number by 


, 2 ^ - 2 . 1/2 
a = <a + B ) 
r r r 


(59) 


The angle between the direction of a^. and the x axis Is given by 


^ » tan ^ (Bj./a^) , 


(60) 


and the phase velocity Is 


c 


ph 


■ w/ct^ 


( 61 ) 


At this point. It Is tempting to form a complex wave number a by (51) 
with the real part given by (59) and the imaginary part by a similar equation 
in terms of and B^> However, this procedure Is valid only when 


- 6r'“r 


( 62 ) 
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and it Is not possible to make this assumption and have the spatial theory 
produce reasonable results. Instead it is necessary to separate the wave 
amplification from its orientation and introduce the new quantity 


. 2 . „2,l/2 

(Oj + Bj) 


(63) 


which makes an angle 


ijl “ tan ^ (64) 

with the X axis. If x is the coordinate in the direction of a and x is the 

r 

coordinate in the direction of a,, we can rewrite (58) as 

q'(x,y,z,t) = q(y) exp (- ci^x) exp [ i(o^x - wt)l . (65) 

It follows that the spatial amplification rate is 

(1/A) (dA/dx) ” - . (66) 

To be more precise, (66) gives the maximum spatial amplification rate for the 
particular choice of there are lesser amplification rates in other 
directions, and of course is itself a function of ij). We see that for 
three-dimensional waves the spatial theory has a difficulty not present in the 
temporal theory: in addition to the wave orientation angle t(), the maximum 

amplification direction (p must be specified before any calculations can be made. 
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The three cases which correspond to (55) in the temporal theory axe: 


> 0 

damped disturbances. 


■ 0 

neutral disturbances, 

(67) 

< 0 

amplified disturbances. 



3. Relation Between Temporal and Spatial Theories 

We now have a temporal theory In which the computation of the eigen- 
values Is straightforward, but which does not yield a spatial amplification 
rate, and a spatial theory which yields a spatial amplification rate, but only 
after the unknown angle has been specified. The problem of choosing is 
avoided only in the special case of a two-dimensional disturbance in a two- 
dimensional boundary layer where both ilf and i(; are zero. The resolution of 
both of these dilemmas is provided by introducing the powerful concept of 
group velocity. 

A laminar boundary layer is a dispersive medium for the propagation of 
instability waves. That is, differtr^ frequencies propagate with different 
phase velocities, so that the individual harmonic components in a group of 
waves at one time will be dispersed (displaced) from each other at some later 
time* An overall quantity, such as the energy density or amplitude, does not 
propagate with the phase velocity, but with the group velocity. Furthermore, 
the group velocity can be ccnsidered a property of the individual waves, and 
to follow an individual frequency we use the group velocity of that frequency. 
Consequently, an observer travelling with the group velocity of a parti lar 
frequency will always see that frequency and its associated amplitude. These 
concepts were originally developed for fully dispersed wave trains in a homo- 
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geneouB medium with no dissipation (real a, B, oj). However, If the Inhomo- 
geneity of the E*'?diuffl and the wave attenuation or ampliflcatio are both "small" 
over a wavelength, then the concepts still apply. The meaning of "small" can 
be made more precise by a multiple-scale pertu*'bation analysis. These conditions 
appear to be satisfied for Tollmien-Schlichting waves in moderately unstable 
boundary layers. In addition, an initial arbitrary waveform quickly becomes 

the dispersed wave train of the theory because, as shown by Mark^ ' for the 

( 36 ) 

temporal theory and by Comer, Houston and Ross'' ’ for the spatial theory, an 

other modes except tb'> fundamental (Tollmien-Schlichting) mode are heavily 

damped. The basic ideas of linear dispersive wave theory for conservative 

( 37 ) 

systems are thoroughly discussed by Whitham. An application to the non- 

( 38 ) 

conservative boundary layer has been made by Landalil^ ^ (but see also the 

( 39 ) 

criticism by Stewartson ). 

The dispersion relation is 

0) * tt)(a, 6, X, z) , (68) 


and the components of the (vector) ^ •'^up velocity In the x and z directions 
are obtained by differentiating (68) with respect to a and 6- When a, 8 and 
0 ) are real. 


c 

g 


( 3oi Bu) 
3a ’ 38 



(69) 


This same expression can be used in the tempi .al theory with ou replaced by 
and in the spatial theory with ot^ and 8^ oi and 0, The imaginary 
part of the group velocity is neglected (and is zero at the point of maximum 


amplification rate). 


- 19 - 



77-15 


From uhat has already been said, it is clear that the temporal and 
spatial amplification rates are related by the group velocity. That is, from 
the parallel-flow taaforal viewpoint we can form a spatial amplification rate 
by following the wave with the group velocity (now independent of x and s) . 
The time derivative is thus transformed into a space derivative by 


d d 

— • c- — . 

dt * dx 


where x must be in the d..rection of c . Consequently, 

8 


(a>. 


^ |c !" * ‘ 

' g’ 


and the direction of can be written 


= tan ^ I (3u^/3B)/(8ti)^/3a)] . 


(70) 


(71) 


(72) 


The group velocity referred to here is the group velocity of the temporal 
theory . 

The prcblem of converting a temporal to a spatial amplification rate was 

(4) 

first encountered by Schiichting , who used the two-dimensional versions 
of (71) and (69) without comment. The same relation was also used later by 
Lees^^^\ but the first mathematical derivation was given by Caster for 
the two-dimensional case. The derivation of Caster's relation is quite 
simple and can easily be generalized to three dimensions. The derivation 
starts from the fact that the complex frequency u is an analytic function of 
the complex variables a and 6 at a fixed x an' Therefore, the Cauchy- 
Riemann equations 
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( 73 ) 


can be applied. In the two-dimensional case, the left-hand side of the f '"st 
of these equations is the group velocity. The derivative on the right-hand 
side can be approximated by noting Chat decreases from its temporal value 
to zero in the spatial theory as goes from zero to its value in the spatial 
theory. If the amplification rate is small, these variations can be considered 
linear, and 


3 ( 1 ) 




(74) 


We see that (74) is the same as the two-dimensional form of (71) with 
the Important difference that the relation is revealed to be only an approxi- 
mation valid for small a^. In the three-dimensional case. If ip is specified 
arbltiarlly and the x axis rotated to lie in the ip direction, (74) will still 
apply with a. replaced by a. and c by the component of the group velocity 

X X g 

in the (p direction. Consequently, when i is chosen in Che direction of the 
group velocity, (71) immediately follows, but again as an approximation rather 
than as an exact expression. 


C. Numerical Procedures 

1. Types of Ifethods 

(3) 

Since the early 1960's, the asymptotic theory developed by Tollmien 
(4) 

and Lin has been largely superseded as a means of producing numerical results 
in favor of direct solutions of the governing differential equations on a 
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digital computer. The numerical methods employed fall roughly into three 

(42) 

categories: (1) finite difference methods, first employed by Thomas in 

his pioneering numerical work of 1953; (2) shooting methods, first employed 
by Brown^^^^ (a successful low Reynolds nundier program was operating as early 

as 1934); and (3) spectral toethods, first used by Gallagher and Mercer^ ^ 

(44) 

with Chandrasekhar and Reid functions^ and later improved by Orszag with 
the use of Chebyshev polynomials* All of these methods have advantages and 
disadvantages which show up in specialized situations* but all are probably 
equally able to do the routine eigenvalue computations required in transition 
prediction calculations* However* it is the shooting methods that have mainly 
been applied to this problem and will be described here* 

(45) 

After Brown's initial work* programs were developed by Mack* 

Landahl and Kaplan* Radbill and Van Driest Lee and Reynolds *^^^^ 
Wazzan* Okamura and Smith* and Davey*^^®^ among others* Most of these 
programs solve only the Orr-Sommerfeld equation; exceptions are the 

(14) (45) 

compressible program of Brown* and the program of Mack which was also 
originally developed for compressible flow and only later extended to 
incompressible flow* All of the programs except Brown's have the common 
feature that the numerical integration proceeds from the freestream to the 
wall. 

The early applications of shooting methods suffered from the problem of 
parasitic error growth* This growth arises because of the presence of a 
rapidly growing solution (the local "viscous" solution) which any numerical 

roundoff error will follow. The rapidly growing error eventually completely 
contaminates the less rapidly growing solution* Tlie essential advance in 
coping with this problem* which had previously limited numerical solutions to 
moderate Reynolds number* was made by Kaplan. The Kaplan method "purifies” 
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the contaminated solution by filtering out the parasitic error vdienever it 

becomes large enough to destroy linear independence. 

2. Gram-Schmidt Orthonormalisation 

An alternative method, first employed by Bellman and Kalaba^^^^ and 

applied to the stability problem by Radbill and Van Driest and Uazzan, 

(17) 

Okamura and Smith, is that of Gram-Schmidt orthonormalization. This 

method has the advantage that it is easier to generalize to higher-order 

systems thaix is the Kaplan filtering technique. However, the geometrical 

argument often adduced in its support that this procedure preserves linear 

independence by keeping the solution vectors orthogonal cannot be correct 

because the solution vector space does not have a metric. Instead, the method 

works on exactly the same basis as Kaplan filtering: the “small'* solution is 

replaced by a linear combination of the “small" and “large" solutions which 

is itself constrained to be "small." 

For the simplest case of a two-dimensional wave in a two-dimensional 

boundary layer, there are two solutions, and each consisting of four 

conq>onents. In the freestream, Z^^^ is the Invlscid and Z^^^ the viscous 

(3) 

solution. Although this identification is lost in the boundary layer, Z 

continues to grow more rapidly with decreasing y than does Z^^^. The 

( 3 ' 

parasitic error will follow Z ' , and when the difference in the magnitudes of 

Z^^^ and Z^^^ exceeds the computer word length, Z^^^ will no Icnger be 
(3) 

independent of Z . Well before this occurs, the Gram-Schmidt orthonormaliza- 
tion algorithm is applied. The "large" solution is normalized component 

by component to give the new solution 
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triiere the asterisk refers to a complex conjt^ate and {? to a scalar product. 
The scalar product of and is used to form the vector 

s<» . [ ) s<»l /{§<»* s<‘> < 


to replace , uhere S refers to. the quantity in the preceding square 
brackets. 

The numerical integration continues with and in place of Z^^^ 

and Z^ and when in turn |S exceeds the set criterion of« say !0^ with 
single precision arithmetic and a 36 bit computer word, the orthonormalization 
is repeated. With homogeneous boundary conditions at the wall, it makes no 
difference in the determination of the eigenvalues whether the Z*s or S's are 
used. A linear combination of the two solutions satisfies the f(G) ” 0 
boundary condition, but the (^(0) •• 0 boundary condition will in general not 
be satisfied unless a, 3 and <i> satisfy an eigenvalue relation. 

3. Wewton-Raphson Search Procedure 

The Newton-Raphson method has been found to be satisfactory for (Staining 
the eigenvalues. In the spatial theory with oi and 3 fixed, the guess value of 
is perturbed by a small amount (A « 0.001 a^) and the integration repeated. 
Because <|)(0), the third component of the linear combination of the two 
Independent solutions Z^^^ and (or and which satisfies f(0) * 0, 

is an analytic function of the complex variable a, even after orthonormalization, 
tl° Cauchy-Rieraann equations 


3a, 3a * 

i r 

3a^^ 3a^ * 


( 77 ) 
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can be applied to eliminate the need for a second integration with 
perturbed. 

The corrections and 6a^ to the initial guesses and are 
<d>tained from the residual ^(0) and the numerical (linear) approximations to 
the partial derivatives by the equations 


3t_(0) (0) 

■ 5 ^ *“t “5^ *“i ■ - ♦!<“> 


(78) 


Hj(0) H.(0) 

___ 5^^ + - - ♦^(O) 


The corrected and are used to start a new iteration, and the process 
continued until 5a^ and have been reduced below a preset criterion. 

4. A Numerical Example 

As an example we will consider that the Reynolds nund>er and frequency 
are specified for the flat-plate boundary layer, and we wish to determine 
both the complex wave nund>er of the spatial theory and the wave number 
and amplification rate of the temporal theory. It is convenient to define 
the length scale as 


* 

L 


(v X /up 


(79) 


With this choice, the Reynolds nund)er is 


it it . it 

U^ L /v 


* *. *,1/2 

(Uj X /v )^'- 


Re 


1/2 


(80) 


The Reynolds numbers and wave numbers based on the displacement, momentum and 

★ 

boundary-layer thicknesses are obtained by multiplying R and a based on L by 
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the appropriate diieensionlesa thickness. For the flat-plate boundary layer, 
these factors are, in the above order, 1.7208, 0.66411 and 6.0114 (u/Uj^ “ 0.999). 
The usual form of the dimensionless frequency is 

* * 

F » u V /Uj = (i»/R . (81) 

For the example, we choose 

R » 1000, F = 0.3 X lO"^. 

and for the initial guess of the complex a in the spatial theory, 

= 0.12 , » 

and for Che temporal theory, 

a * 0.12 , w. * 10.0 

r 1 

The initial conditions of the two independent solutions are evaluated from 

Che formulas of Section II-A.4, and the numerical integration started at 

* * ^ ^ 
y^^ “ y /L “ 8.0. Each solution is Integrated to the wall (at y ® 0) by 

means of a fourth-order, fixed step size, Runge-Kutta integration. Other 

integration methods are also satisfactory, but this simple method has been 

found to be trustworthy in a wide variety of problems, particularly in the 

(35) 

difficult problem of determining the eigenvalue spectrum. Variable step- 

size Integrators are not recommended. 
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The operation of the eigenvalue search procedure is shown in Table I. 

The search is continued until 6|a|/|aj is reduced below a preset criterion 
(0.005 is usually adequate for three-place accuracy) in the spatial theory, 
and 6o/a and ^ |a>^( / joi^I separately satisfy the criterion in the temporal 
theory* Since < 0 ^ is held constant in the temporal theory, two perturbation 
integrations are required per iteration instead of one as in the spatial theory. 
Two iterations are usually sufficient to achieve convergence in a large scale 
computation where previously calculated eigenvc jes can be used to make good 
Initial guesses. Poor initial guesses were deliberately chosen in the examples. 

Once a single eigenvalue has been found for a given boundary layer, all 
others can be readily obtained. Automatic procedures can be Included in the pro- 
gram to produce all of the unstable eigenvalues of a given frequency or wave 
number mesh up to some specified large Reynolds number in a single computer run 
of a few minutes for a two-dimensional boundary layer. The problem of obtaining 
the initial eigenvalue can prove troublesome. If the boundary layer under 
investigation is not far removed from one for which the eigenvalues are known, 
it is always possible to make a close enough guess for the Newton-Raphson 
procedure to converge. If not, it is necessary to perform individual integra- 
tions in the complex a plane (spatial) or oj plane (temporal). Then the contour 

line method of Mack,^^^^ or the more elaborate methods of Antar^^^^ and Jordinson 
(52) 

and Caster, can be used to locate the initial eigenvalue. 

D. Some Numerical Results 

1. Amplification Properties of Falkner-Skan Profiles 

For a two-dimensional incompressible boundary layer on an impermeable 
surface of zero curvature and with no body force or surface heating, the only 
flow parameter left is the pressure gradient. The effect on the amplification 
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properties of systematic changes in the pressure gradient can best be 
demonstrated with the Falkner-Skan family of similar velocity profiles. 

The single parameter of a Falkner-Skan boundary layer is the dimensionless 
pressure gradient 


m 


* * 
X dp. 


2 

* *^ * 
p dx 



It is usual to replace m by 


(82) 


3 = 2m/(m + 1) (83) 

(not to be confused with the lateral wave nundier). Ihe values of 8 range from 
1.0 for the two-dimensional stagnation point boundary layer, through 6=0 for 
the flat plate, to 6 = 0.1988377 for the separation profile. More negative 
values of 6 represent reverse-flow profiles. 

Two theoretical principles are useful in interpreting the numerical 
results. The first comes from the viscous asymptotic theory, and states 
that increasing the negative curvature of the velocity profile near the wall 
Increases the stability. The second comes from the inviscid theory of Rayleigh 
and states that as the inflection point moves away from the wall, the profile 
becomes more unstable. We need only recall th.it as 6 increases from zero, the 
first effect occurs; as it decreases from zero, the second affect occurs. 

Some properties of Falkner-Skan profiles are summarized in Table II, 

ic 

wher^ are listed 5,6 and 5 (the displacement, momentum and boundary-layer 

* ic 

thicknesses) made dimensionless with respect to L ; H = 5 /O, the shape 
factor; U"(0) (■ - m) , the curvature at the wall; Fg/fi* the location of the 
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Inflection point; and c , the dimensionless velocity at the inflection point 

8 

and equal to the phase speed of the neutral inviscid disturbance. 

The temporal stability of Falkner-Skan profiles was computed in great 
/o\ 

detail by Pretsch^ ^ from the asymptotic theory (these results may also be 

(53) 

found in Smith and Gamberoni )» and both the temporal and spatial stability 

(17 18) 

from direct numerical solutions by Wazsan, Okamura and Smith. * For 
application to the transition problem, the quantity of primary interest is 
the integral of the spatial amplification rate for a constant frequency. With 
L still defined by (79), R by (80) and by (82), it follows from (66) that 



o 

where is the amplitude at the initial Reynolds number R^. It is convenient 
to take R^ as the lower-branch t atral stability point (initial point of 
instability for the frequency under consideration) In order to give a unique 
meaning to A/A^. If the definition of the dimensionless frequency given by 
(81) is to be retained, then F is no longer constant for a constant dimensional 
frequency, but is given by 


F(R) 


it it 
(A) V 


* 

Ui (x ) 


it it 
(a V 


* / 

"l <’'c > 



4m/m+l 


(85) 


Thus the calculation starts at = U, (x ) x /v with dimensionless frequency 
- o 1 o o 

it it it^ it 

F * 01 V /U- (x ) , and the integral (84) is evaluated using a.(R) calculated 
o 1 o 1 


with the F(R) from (85). 
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‘Hte length scale L is not the only possible choice with which to make 

fg ^ if 

dimensionless. The inverse of the unit Reynolds number, v /U^, is often 
found in the literature. With this choice. 


In 




nH*l 



d Re 


( 86 ) 


o 

where Re is the x-Reynolds number. This procedure is perfectly acceptable, but 
it has the disadvantage that a^v /\J^ (= a^/R) always goes to zero as Re ® 
even when the boundary layer is unstable to inviscid disturbances. On the 
contrary, the quantity is based on a boundary-layer length scale, and so 
can also be used as the inviscid amplification rate. Comparisons are easy 
to malce between the viscous and inviscid stability theories in terms of 
a.. For example, with 3 “ - 0.15, the inviscid (-a.) is 0.0199, and in 
the viscous theory already reaches the slightly higher value of 

0.0202 at the low Reynolds number of R « 400. 

For the transition-prediction calculations of Section IV, three 
quantities will be needed: (i) the envelope curve. In (A/A^)^^ vs. R, 

formed by the individual In (A/A^) vs. R curves; (ii) the frequency, » 

the envelope curve; and (3) a bandwidth AF of the frequency response 
curves to be discussed below. The envelope curves for sever'^l values of 3 
are shown in Figure 1. Both the strong stabilizing effect of a favorable 
pressure gradient, and the even stronger destabilizing effect of an adverse 
pressure gradient are clearly evident. The corresponding frequencies are 
shown in Figure 2. (The 3 « - 0.05 curve is almost coincident with the 
3 “ 0 curve and so cannot be shown.) We see a clear distinction between 
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viscous and inflectional Instability. Viscous instability is primarily a 
low frequency instability, particularly when the instability only develops 
at high Reynolds numbers. In contrast, inflectional instability is a high 
frequency instability. The ratio of the frequency which gives In (A/A^) •» 9 
for 6 = -0.15 to the corresponding frequency for $ ■ 0.20 is 36. 

The need for quantity (ill) is shown in Figure 3, where In (A/A^) for 
the flat-plate b( undary layer is plotted against F for several Reynolds 
numbers. As R increases, the maximum amplitude ratio increases and the band- 
width of \mstable frequencies decreases. This sharpening of the boundary- 
layer response must be taken into accovmt in any transition prediction method 
that attempts to calculate the integrated (over frequency) amplitude in the 
boundary layer. A quantitative measure of this effect is provided by the 
bandwidth 


where the second term is the frequency at which A/A is 1/e of (A/A ) on 

o o max 

the low frequency side of F ^^ ^. It is desirable to define AF this way 

rather than as a two-sided bandwidth because the practical requirements of 

a large scale eigenvalue calculation are such that A/A^ can always be 

obtained for frequencies smaller, not larger, than Figure 4 gives 

the ratio AF/F as a function of R fur the same $ as in the previous 
max 

figures. 

2. Comparison of Temporal and Spatial Amplification Rates 

The requirement to know the direction of the group velocity before 
computing eigenvalues from the spatial theory for other than two-dimensional 
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waves in two-dimensional boundary layers makes it worthwhile to examine how 
accurate are the spatial amplification rates obtained from the temporal 
theory and Caster's relation in this more restricted case. The group velocity 
of the temporal theory, as given by (69), can be computed by numerical 
differentiation along with the other stability properties provided is a 
smoothly varying function of the wave number. 

Table III presents both and for four (three impllfled, including 
the maximum, and one damped) at three different $ and at Reynolds numbers with 
large amplification rates. The case 6=0 provides amplification rates typical 
of boimdary layers with small favorable pressure gradients; 6 = - 0.10 
provides about the largest amplification rates that can be expected in practice; 
and the separation profile caa be viewed ar providing an upper limit for 
boundary-layer amplification rates, but not particularly representative of 
actual practice. (The maximum amplification rate of this profile occurs at 
R ->• oo and is -0.0480, which is only 12% larger than the R = 300 value of 
-0.042A.) 


From Table III we see that the Caster relation is satisfied quite well 
at 6=0, and less well at the two other 6*s. At 6 = 0, the frequency and 
phase velocity from the two theories also agree closely. It can also be 
noted that the ratio - a>./a. is less than the group velocity at the maximum 

X X 

ampi i.i icatioc rate for all three values of 6. The maximum difference between 

Cg and for amplified disturbances occurs at the maximum amplification 

rate, and is 1.6% for 6=0, 2.8% for 6 “ -0.10, and 5.7% for the separation 

profile. In the latter case, there are also Important differences in and 

c . between the two theories at both low and high wave numbers. The conclusion 
ph 

can be drawn from Table III that the temporal theory and Caster's relation 
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of'er a satisfactory method of obtaining spat:*al amplification rates for 

zero and favorable pressure gradient s, but that this approach becomes 

increasingly more unreliable as the adverse pressure gradient increases. 

3. Group Velocity Direction for Oblique Waves 

The importance of the choice of the direction of maximum amplification 

in the spatial theory can be readily demonstrated by means of oblique v/aves 

in the flat-plate boundary layer. Table IV gives results for the two- 

-4 

dimensional wave of dimensionless frequency F»0.3xl0 atR= 1600, 
and for three oblique waves with rp « 45^, 60^ and 75^. The Reynolds numbers 
and F for these waves were chosen so that in the ^ direction the Reynolds 
number (ft * Rcos^Ji) and dimensionless frequency (F « F/cos i|>) are the same as 
for the two-dimensional wave. According to (22), the Squire transformation, 
the spatial amplification rate in the x-dlrection is 

^i * ^i ^ ’ 

where is the spatial amplification rate of the two-dimensional wave with 

-4 ^ -3 

F ® 0.30 X 10 at R ^ 1600. Since » -3.82 x 10 from Table IV, the 

other three values should be -2.70, -1.91, and -0.989 for = 45^, 60°, 75° 

respectively. We can see that, as already stated in Section II-A.2, this 

relation is true only when For if < is larger than these values 

by amounts ranging up to 39%. Consequently, it is not ^rmissible to use 

^ in order to preserve the real form of the Squire transformation. 

The correct ijl is the direction of the group velocity. Figure 5 

st^s how this direction»now identi||,ed as \{; and computed from the temporal 

theory, varies with for two fixed values of o at R * 1600. The important 
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conclusion to be drawn from Figure 5 is chat regardless of the orientation of 
the constant phase lines (the crests), the group lines, which give the 
direction of energy propagation, remain concentrated near the freestream 
direction. Furthermore, the 10^ entry of Table IV sh''ws that is nearly 
the same as with i|) » 0. Thus, for oblique waves in a two-dimensional boundary 
layer, it Is a permissible approximation to use Che spatial theory with ii « 0^. 
However, it is still necessary to use either the complex form of the Squire 
transformation or the formulation given in Section I1-A.4 to compute the 
eigenvalues . 


III. STABILITT EXPERIMENTS 
A. Schubauer-Skraastad Experiment 
The linear stability theory long went unapprer iated except by its 
founders because of the lack of any convincing experimental confirmation. 

This needed confxrmation was brilliantly supplied by the new classic experi- 
ment of Schubauer and Skramstad^^^ which vas carried out at the National 
Bureau of Standards in the early 1943* s, but not published because of war- 
time restrictions until 1948. This experiment owed its success both to the 
ingenuity of Che experimenters and to the development for the first time of a 
wind tunnel with a really low turbulence level (about 0.03% in the working 
section). When, with this low turbulence level, the signal of a hot-wire 
anemometer placed in the boundary layer of a flat place was displayed on an 
oscilloscope screen, modulated sinusoidal wave trains with almost no random 
character were clearly seen. Schubauer and Skramstad demonstrated conclusively 
that these were true boundary-layer oscillations, and that they were the 
cause, and not the effect, oi transition. 
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In order to make a more quantitative connection with the theory of 

Tol?mien and Schlichting, they used a vibrating ribbon in order to produce 

disturbances of a fixed frequency with a controlled initial amplitude- The 

hot wire then measured the wave length « phase velocity, and amplitude of 

the artificially produced waves as a function of Reynolds number. Numerous 

comparisons were made with the theory which were on the whole satisfactory, 

although the asymptotic theory at that time did not yield very accurate 

numerical results* However, it is well to emphasize that the virtue of this 

experiment did not lie in an exact quantitative correspondence with theory, 

but rather in the systematic way that all essential features of the theory 

were shown to be correct, and in the way the oscillations were shown to 

be necessary precursors of transition- 

The hot wire measures directly the rms disturbance amplitude as a 

function of downstream distance, and the amplification rates must be deduced 

from the slopes of such measurements. Although the interpretation of the 

measured amplitude has difficulties associated with non-parallel flow 

effects, it is still of considerable interest to compare the amplitude 

measurements with the quantity A of stability theory. The comparison is 

given in Figure 6. The experimental wave amplitudes (u* at a fixed distance 

from the wall) are all referred to the amplitude at = 2 in. behind the 

ribbon which was located 4 it. behind the leading edge of the plate with 

Uj * 64 ft/sec. Therefore, the integration of was also started at x^ 

(R “ 1256) , and A is the amplitude at x rather than at the neutral 
o o o 

stability point. The experimental points in Figure 6 show conclusively that 
the frequency is the fundamental parameter that determines whether a wave 
will be amplified or damped, and the agreement with theory is satisfactory 
although far from exact. 
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This is a good place to bring out some import a* quantitative aspects of 
instability waves. The ratio of the wavelength to t..-*. displacement thickness 
is 

T*-|Vh72^ • <“> 

6 \ 1 « Re^ I 

At R 1256 (the of Figure 6), * 0.128 for the frequency with the 

largest amplification rate (f = 120 hz, F = 0.311 x 10*^). Hence \ * 28.5 6* 

(« 8.2 6), and we see that unstable Tollmien-Schlichting waves are long 

compared to the boundary-layer thickness. This result from a theory that 

assumes parallel flow was one of the original criticisms made of the theory, 

'nd also explains why a considerable effort has been made recently to develop 

(32 33,54) 

non-parallel theories. ’ ’ Another interesting point is rbp rapidity with 

which the waves grow. The 120 hz wave increases in amplitude 2.5 times in 
* 

6 in. (» 88 6 , or 3.1 X). From the definition of the amplification rate, 
the fractional change in amplitude per wavelength, with and assumed 
constant, is 


6A 

“T * “a ’ 
r 


(89) 


or about 33% for the 120 hz wave. Although this large growth only exists 
near the maximum amplification rate, it still raises problems concerning 
the application of kinematic wave concepts to instability waves. 

B. Other Older Experiments 

The next stability experiments were carried out by Liepmann.^^^'^^^ 
These experiments were designed primarily to study the stability and 
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transition of boundary layers on concave surfaces » i.e. Gortler instability* 
but measureiaents were also made on the convex side of the curved plate. Two 
plates were used* with radii of curvature of 20 ft. and 2-1/2 ft.* but only 
the 20 ft* plate was used for the detailed measurements on the convex side. 
The boundary- layer velocity profile was close to that of a flat plate* but 
the neutral-stability curve was found to define an unstable region somewhat 
larger than in the Schubauer-Skramstad experiment* perhaps reflecting a 
slight average adverse pressure gradient. Transition measurements on the 
convex sides of both plates showed no effect of curvature* but the imposition 
of favorable and adverse pressure gradients on the 20 ft. plate did produce 
an effect* The transition Reynolds ntimber was increased for a favorable and 
decreased for an adverse pressure gradient* and the percentage change was 
greater for the favorable pre ^ gradient* No stability measurements were 
made with a pressure gradier«^ indeed nc theoretical stability results 

on the effect of pressure gra *ere available to Liepmann at the time of 

the experiments* 

The next experiment* by Bennett studied the influence of freest ream 

turbulence on the instability of the flat-opiate boundary layer. A grid was 

installed upstream of the test section to raise the turbulence level to 

0.42%. The signal from a hot wire in the laminar boundary layer shoved 

fluctuations of a random nature resembling turbulence with little evidence of 

C58) 

Tollmien-Schlichting waves* an observation in accord with Dryden*s work 
of 20 years previously* However* when the power spectrum of the fluctuations 
was measured with a wave analyzer* a local peak was found to develop at 
about the most unstable frequency of linear stability theory* This peak 
grew with increasing dotmstream distance* and gradually disappeared after 
transition started as the spectrum evolved into one typical of a turbulent 
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boundary layer* Measuresients were made at too few x stations in the laminar 
region to determine the amplitude history of the peak« and consequently no 
conclusion is possible from this experiment on the crucial point of whether 
the linear disturbance growth is affected by the external turbulence level* 

A flow visualization technique for instability waves using tellurium 

( 59 ) 

coated rods was developed by Wortmann for use in a water boundary layer. 
When a voltage was applied to the rods, tellurium ions, which are black, were 
released and made the instability waves clearly visible. It was then possible 
to make quantitative measurements, and a neutral stability curve was determined 
that compared favorably with Schubauer and Skramstad*s. Other flow visualiza- 
tion techniques that have been used are smoke, dye^^^^ and hydrogen 
( 63 ) 

bubbles.^ ^ These methods all show the presence of Tollmien-Schlichting 

waves, but have been applied mainly to study the details of the transition 

process following linear amplification. 

The hot-wire anemometer continued to be the primary tool of the NBS 

transition studies, which, following the work of Schubauer and Skramstad, also 

concentrated on the non-linear region. Schubauer and Klebanoff studied 

the characteristics of the turbulent spot; Klebanoff and Tidstrom^^^^ and 

( 66 ) 

Klebanoff, Tidstrom and Sargent^ ^ the sequence of events from the end of 
the linear region to the first appearance of a turbulent spot. It was found 
that the initially two-dimensional Tollmien-Schlichting wave develops a span- 
wise periodicity in amplitude xdiile it is still undergoing linear amplifica- 
tion. There appears to be a characteristic spanwise wavelength that has 
since been observed in other wind tunnels, but no convincing explanation has 
yet been given for the origin of this periodicity. The Knapp and Roache^^^^ 
smoke pictures clearly show that a similar three •dimensionality also 
develops in natural transition. 
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C. Three Recent Experituents 

1. Ross> 3ames> Bnras and Ross Experiment 

What azoants to a repetition of that part of the Schubauer-Skramstad 

experiment which measured the instability of a flat-plate boundary layer was 

carried out by Ross> Barnes* Bums and Ross.^^^^ Their wind tunnel had the 

same turbulence level as the NBS tunnel, was of sisiler size, the flat plate 

was also mounted vertically, aud a vibrating ribbon was used to produce the 

instability waves. The hot-wire measurements were compared with the detailed 

f 681 

calculations of Jordinson,^ ' and excellent agreement was obtained for the 

distribution of u* through the boundary layer. Quite good agreement was also 

obtained for the A/A^ of three frequencies (F x 10^ » 0.82, 1*10, 1.57'. 

Particular attention was paid to the region of the minimum critical Reynolds 

number and the maximum unstable frequency. Measurements in this region are 

difficult because the boundary layer is both thin and rapidly growing. A 

—A 

neutral-stability curve was arrived at with R = 230 and F - A x 10 
^ cr p 

A 

compared to the theoretical values of R = 302 am F * 2.5 x 10 . The 

cr p 

differences between theory and experiment were attributed to non-parallel flow 

effects, a supposition since confinmid by the good agreement of the Saric- 
(33) 

Nayfeh theory with the above experimental results. 

It must be pointed out that the frequencies that define the 
portion of the neutral-stability curve are not those that are important to 
transition in environments with small disturbances. For example, at the 
Schubauer-Skramstad transition Reynolds number of 2.8 x 10^, we see fron 

•4 

fiaure 2 that F * 0.29 x 10 compared to F = 3.5 x 10 at R = 2r>0 
according to Saric and Mayfeh. Tlie former frequency first becomes un?itable 
at R * 820, where the non-parallel flow effects are less severe than at the 
minimum critical Reynolds number. The non-parallel flow effects will be 
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even further reduced in boundary layers with favorable pressure gradients 
where instability occurs at much higher Reynolds numbers than for the flat 
plate. 

2. Klebanoff-Tidstrom Experiment 

It is surprising that tne traditional low-speed stability theory has 
never been tested eKperimentally for other than the Blasius boundary layer 
studied by Schubauer and Skramstad. There have been boundary- layer stability 
experiments on the effects of rotation, compressibility, free 
convection, and a heated wall in water, but the accompanying 
theories represented extensions of the existing theory to new flow situations* 
There have been transition experiments on the technically important effects 
of pressure gradient and suction, but no stability experiments. 

There is an important aspect of transition that also has received 

little attention, and that is the relation of a particular disturbance 

source to the transition process in a boundary layer. In other words, the 

precise mechanism by which, say, freestream turbulence, sound, and different 

types of roughness cause transition remains to be discovered. Only in the 

case of two-dimensional roughness has the mechanism been found thanks to a 

(73) 

remarkable experiment by anoff and Tidstrom. For circular roughness 

ie 

elements with a diameter of about 0.8 6 , these investigators determined that 
transition was moved forward from its normal location, not by a disturbance 
introduced into the boundary layer by the roughness, but by the increased 
amplification of an already existing instability wave in the pressure 
recovery zone behind the roughness. In a certain sense, w3 do have here 
an example of a stability experiment with a pressure gradient, but this 
very special flow is scarcely representative of pressure-gradient boundary 
layers. 
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Although the demonstration by Klebanoff and Tldstrom of the actual 

transition mechanism was complete in itself and did not rely on theoretical 

comparisons, they did compare their disturbance growth measurements with 

( 8 ) 

results obtained from Pretsch's charts. ^ Unseparated adverse pressure 
gradient Falkner-Skan profiles were fitted to the measured profiles on the 
basis of the shape factor H. The calculated and measured growths were in 
good agreement. However, there are three objections that can be raised 
to Klebanoff and Tidstrom*s procedure; (i) the profiles in the Initial 
part of the recovery zone are clearly separated; (ii) the Pretsch charts 
are not very accurate for adverse pressure gradients; and (iii) the group, 
and not the phase, velocity should have been used to transform temporal to 
spatial amplification lates. 

When the computation is repeated with the correct spatial amplification 
rates for the fitted Falkner-Skan profiles, the results are not in agreement 
with the measurements. The obvious next step is to use the correct velocity 
profiles which were measured in great detail and with little scatter. 
Unfortunately, the separated-flow region, which has a strong influence on 
the instability, could not be measured because of its closeness to the wall. 

In any case, stability calculations based on experimental curve fits to 
profiles with an inflection point are not likely to be meaningful. The 
conclusion to be drawn is that although the Klebanof f-Tldstrom experiment 
firmly established that two-dimensional roughness influences transition by 
destabilizing the boundary layer, the amplification measurements cannot be 
used as a test of a particular form of the stability theory in a rapidly 
varying pressure-gradient flow. 
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3. Gaster^Grant Experiment 


Another recent stability experiment of great interest was carried out 

(74) 

by Caster and Grant In the same wind tunnel used In the Ross» Barnes » 

Bums and Ross e^erlment.^^^^ All of the other stability experiments with 
artlficJally produced disturbances have followed directly in the Schubauer- 
Skramstad tradition of dealing with a Tollmien-Schlichting wave of a single 
frequency. This approach was followed even where the traditional vibrating 
ribbon was replaced by a siren^^^^ or a glow-discharge device. In the 

latter instance» oblique waves were also produced^ but still only of a single 
frequency. Caster and Grant used the completely different approach of trans- 
mitting an acoustic pulse from a loudspeaker through a small hole in their 
flat plate to produce a pulse disturbance in the boundary layer. This 
disturbance, after the rapid decay of all higher modes consists of Tollmien- 
Schlichting waves of all frequencies and orientations. This type of experi- 
ment, while not as suited to mapping out stability boundaries and making 
the other usual checks of stability theory as in a vibrating ribbon experiment 
is in some respects closer to the transition problem. An external disturbance 
such as freestream turbulence or sound, will also produce Tollmien-Schlichting 
waves of all frequencies and orientations, but in a random manner that makes 
the ensuing wave motion in the boundary layer difficult to sort out. The 
pulse experiment gives a controlled disturbance of this type which enables 
the details of a group of waves (a wave packet) , rather than a single 
Fourier component, to be investigated. 

In the experiment, the hot wire was placed at 13 spanwise stations 
at each of six downstream locations to record the passage of ^he wave 
packet. Although several theoretical treatments of the motion of a pulse 
in a shear layer have appeared, these papers all made use of asymptotic 
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methods. In order to have more exact numerical results, the amplitude of 
the wave packet was calculated from 


A(x,z,t) 


o> 

II 


-1/4 ( 

X exp ; 


1 / 


o(x) dx + 6z-a)t I dB db) 


(90) 


where x^ Is the location of the pulse; and the Fourier coefficient has been 

-1/4 

set equal to uni-.y. The factor x Is Intended to account for the 
effect of the boundary-layer growth, end the Integral of a(x) appears 
for this same reason. The eigenvalues a(x) were obtained from the 
spatial stability theory for all values of B and ta by use of the Squire 
transformation and a complex Reynolds nundier. In order to evaluate the 
double Integral with sufficient accuracy to reproduce the wave packet, 
10,000 eigenvalues were needed. The agreement that was achieved 
between theory and experiment for about the first 2/3 of the distance 
covered was Impressive. After this, the experimental wave packet 
distorted In a way nor given by the theory, possibly due to non-llnear 
effects. 
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IV. TRANSITION PREDICTION 
A. Nature of the Problem 

A boundary layer hap a specific displacement thickness and skin friction, 
but it does not have a specific transition Reynolds number. The observed transi- 
tion Reynolds number depends on the presence of disturbances in the boundary 
layer, which in turn are related to various disturbance sources. If there 
were no disturbances, there would be no transition and the boundary layer would 
remain laminar. Consequently, it is futile to talk about transitioii without 
in some way bringing in the disturbances which cause it, and any transition 
criterion of an empirical nature can only be valid for a very specific dis- 
turbance environment* 

Another point which must be settled before going deeper into the subject 
is to define the circumstances under which linear theory can be used to predict 
transition. Disturbances such as a large three-dimensional roughness element, 
or an air jet, cause transition to occur in the immediate vicinity of their 
location. However, In many other instances, the disturbances act in a more 
indirect manner. A moderate freestream turbulence level or acoustic intensity 
affects transition by producing Tollmien-Schlichting waves in the boundary layer 
which then amplify, distort, and finally culminate in the sudden appearance of 
a turbulent spot.^^^^ These instability waves can be described initially by 
linear theory, and if the initial disturbance amplitude is sufficiently small, 
the region of linear growth can be of significant extent. In this case, the 
exponential growth of a linear disturbance, and the absence of an extensive 
region of dis*:orted laminar flow before the sudden breakdown to turbulence 
that is such a distinguishing feature of low-speed boundary- layer transition, 
makes it feasible to base a method of transition prediction directly on the 
linear theory itself. 
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B. Anq>litude Density Methods 


1. Liepmann*s Method 

(U) 

The first application of linear theory* to transition, by Schllchtlng, 
has already been mentioned In Section I. A later method, developed by 
Schllchtlng for application L'> airfoils and described In his book,^^^ was 
based on the minimum critical Reynolds number, and thus avoided having to 
Introduce the disturbance level. It was Llepsiann^^^^ who first worked out 
a formula that Included most of the ingredients needed to use linear theory 
in transition prediction. This method can be classified as an aisplltude 
density method, because It considers only a single frequency component, i.e., 
a single spectral line, from what must be a continuous disturbance power 
spectrum. 

Llepmann's idea was that transition should occur when the ratio of 
the Reynolds stress of the boundary- layer disturbance equals the mean viscous 
stress, and that the linear theory provides an adequate means of computing the 
Reynolds stress. Llepmann's formula is 


T_ 

T 


C, 


L Jo' 'o 


(91) 


where is the skin- friction coefficient. Liepmann used results of 
Schlichting^^*^^^ to evaluate all terms in (91) except the initial 
amplitude He believed to be related to the freestream turbulence 

level, but the details of this relationship are as unknown today as they 
were 30 years ago. 
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2. Method 

(53) 

The next step was the independent development by Smith and Gainb«ironi 
and by Van Ingen^®^^ of what is often called the e^ method. Smith and Gamberoni, 
starting from (91) , found thac the only quantity on the right-hand side that 
they could compute other than the skin friction, namely the amplitude ratio 
A/A^, was itself sufficient to correlate a large number of experimental mea- 
surements of transition. With the use of the stability charts of Pretsch, 

/ QOX 

and an approximate method of computing the boundary layer developed by Smith, ^ 

transition was found to occur when jEn(A/A ) ^ 9. A similar result was obtained 

o 

by Van Ingen, but with the exponential factor equal to 7 or 8. This method has 

been justly criticized for basing transition on a ratio and not on the dis- 

tutbance amplitude. However, this criticism ignores the fact that the wind- 

tunnel data which were used to develop the method probably refer to rather 

similar disturbance environments, and a'' so doesn't give credit to the value 
9 

of the e method in comparative studies. For a fixed value of A^, which is 
equivalent to a fixed disturbance environment, the disturbance amplitude in 
the boundary layer does in large measure vary as A/A^. Consequently, the 
effect on transition of changing a parameter which governs the mean boundary 
layer, such as the p^^ssure gradient, can in this particular circumstance be 
estimated by means of the single factor 

In more recent work, Jaffe, Okamura, and Smith^ ^ replaced the original 
approximate numerical methods of Smith and Gamberoni with more exact methods , 
and Van Ingen^ ' widened the range of applications. However, the method 
remains essentially as originally developed, and the key to success still lies 
in a judicious choice of the value of the exponential factor. The use of the 
method is simplicity itself once the necessary numerical tools have been 
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assembled. Three computer programs are needed. The first computes the 
inviscld preasure distribution over a specified planar or axis 3 onmetric 
body shape; the second uses this pressure distribution to compute the laminar 
boundary layer; and the third computes the spatial amplification rate and 
its integral, the amplitude ratio. Transition is considered to occur whenever 
jKn(A/A^) reaches the chosen numerical factor. 

The ^ly aspect of the method that remains to be mentioned is the 

equation for A/A^. With a non^similar boundary layer, there is no simple 

relation between the arc length x and the boundary- layer thickness as ^ith 

the Falkner-Skan profiles. One convenient expression is 

* 


in 



•k 

U 


00 


•k 

V 



(92’> 


o 

where c is the chord or body length, is the free.tream velocity, (x ; 

is the edge velocity, is the local displacement thickness Re>nolds number, 

and (of.).* ® 01* 6*. As before, the integrcvl is evaluated for a constant 
1 0 i 

dimensional frequency. Another expression, which preserves the boundary- 
layer length scale L*(x*) = [v* ]^, is 


in 




(93) 


where R * (IJ * x*/v*)^ and ct = (V 
00 00 11 

3. Modified Method 

9 

The main problem with the e method lies in the proper selection of 
the exponential factor, which is by no means always equal to nine. One ^.ay 
to avoid an arbitrary choice is to relate the factor, which can be called n, 


I 
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to the disturoance level in an eioplrical manner. This procedure has been 

applied to the case of freestream turbulence by Van Ingen^^^ and Mack.^^^^ 

WiiKt-tiinnel data on the influence of freestream turbulence on transition on 

f 28) 

a flat pleta have been collected by Dryden^ ^ and are shown in figure 7. 

A sore recent and more complete collection of data has been assembled by 

y oOx 

Hall and 6ibbing8> * but the additional data still follow the trend of 
figure 7. The start-of-traawltion Reynolds mind>er is Re^ and the turbulence 
level is Uj^/Oj^, or T « + w^)/3 If all three fluctuating 

velocity components have been measured. The turbulence level in a low tur- 
bulence wind tunnel is Increased by successively removing the damping screens > 
and high levels are achieved by Installing grids just upstream of the test 
section. The total disturbance level in a wind cunnel is made up of both 
turbulence and sound. Below T « 0.1%, the sound component controls 
transition, its a result, decreasing the turbulence component only decreases 
the s:^gnal registered by a hot-wire anemometer withcmt affecting the transition 
Reynolds number. It is for this reason that the curve in figure 7 is level 
for T <0.1%. The effect of freestream turbulence on transition is given 
by the sloping portion of the curve for T > 0.1%. 

An analyc^ls of results computed from the stability theory shows that 
this portion of the curve can be accounted for by considering A to remain 
fixed at trar..<<ition, and with 

A . (94) 

o 

This same variation in A/A^ Is given by letting the exponential factor vary 
according to 


n * - 8.43 - 2.4 An T . 


(95) 
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This method, where n is related directly to the disturbance level » may be 

9 

called the modified e method. Whether it is valid for other than flat* 
plate boundary layers cannot be determined until systematic transition data 
comparable to figure 7 become available for more gei^ral boundary layers. 

C. Amplitude Method 

1. Formulation of Amplitude Relation 

Although the quantity A has been referred to in the foregoing as the 

amplitude, it is more properly called an amplitude density because we have 

been dealing solely with disturbances of a single frequency. That is, A 

represents only a single spectral line of a continuous power spectrum. The 

9 9 

method of Liepmann, the e and modified e methods are all amplitude density 

methods. An amplitude density method which considers A to be a function of 

the energy density of the one -dimensional power spectrum of the external 

f 89) 

disturbance was developed by Mack^ " for transition in a supersonic wind 
tunnel* Although this method was reasonably successful in explaining the 
effects of Mach number, unit Reynolds number and tunnel size on transition, 
it can be criticized on the basis that a single spectral line is not an 
adequate basis for relating the disturbance spectrum in the boundary layer 
to the spectrum of the external disturbance source. 

Consequently, another method has been developed which approximates 
the actual wide-band disturbance amplitude in the boundary layer* The starting 
point is the following expression for the disturbance amplitude A^ as a function 
of the Reynolds number R = (U^x/v)^: 
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In order to clarify the derivations, the use of asterisks has been dropped 
and all quantities are new dimensional except A, and A^ which are referred 
to the freestream velocity. In (96), A is the length scale of the external 
disturbance, w its (circular) frequency, and 3 its lateral wave number. The 
initial disturbance in the boundary layer is considered to be produced directly 
by the external disturbance and thus must be scaled in terms of A, while the 
amplitude ratio continues to be scaled in tersm of the boundary- la3rer length 
scale. The velocity scale for both A^ ami A/A^ is - ^^e basic assumption 

in (96) is that all phase relations are random so that the various harmonic 
components add in the square. Another viewpoint is that A^ is the long time 
average of the disturbance amplitude, and we will be determining the transition 
Reynolds number of this average disturbance. What we would really like to Know 
is the average transition Reynolds number produced by a disturbance source that 
is steady only in its time average, and it is by no means certain that the two 
Reynolds numbers are the same. 

It is entirely possible to evaluate the double integral of (96) numeri- 
cally once A^ is known, but it is more prudent to adopt a simpler approach and 

9 

keep the numerical requirements nearly the same as for the e method. The two 
limiting cases which can be considered are where the bandwidth of the boundary- 
layer response is small compared to that of the external disturbance, and where 
the opposite prevails. We will treat only the first of these. Hence A/A^ 
acts as a 6- function and (96) can be written as 
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We may further take advantage of che fact tliat the xeaponac curves A/A^ 
are cuten rather siftiilar in shap^, and evaluate the double integral in an 
approximate manner to arrive at 



The last two factors are consistently defined bandwidths of the frequency 

2 

and lateral wavenumber response curves . The factors (U^A/v) and 1/R come 
from converting the integration variables (nA/U^^ and BA into the boundary- 
layer variables (ov/U^^ and The constant expresses the difference 

between the exact and approximate integrations. 

2. Interaction Relation 

In order to proceed further in the evaluation of it is necessary 
to relate A^ to the particular external disturbance under consideration. 
Unfortunately, at the present time nothing is known about the mechanism by 
which any external disturbance produces Tollmien-Schlichting waves. Although 
we could proceed on a purely empirical basis, a better appreciation of the 
problem is achieved by adopting a particular viewpoint. The viewpoint we adopt 
is that the instability waves are produced in the viscous sublayer, or Stokes 
layer, set up at the wall by the freestream disturbances. The forced response 

in this layer to a sinusoidal disturbance is obtained from the simplified model 

(2) (90) 

of Prandtl. This approach was used by Sternberg in studying the viscous 

sublayer of a turbulent boundary layer. The final step is to consider to 

be proportional to the induced normal velocity of the forced response. The 
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idea of relating the free to the forced response had some success in studying 

(91) 

supersonic boundary layers irradiated by strong acoustic disturbances. 

A straightforward analysis leads to 



where vj is the rms nonoal 'Viscous** velocity of the harmonic ccmtponent 

((oA/^t > 3A) > p^ is the corresponding rms pressure fluctuation at the wall, 
i (0 

and t is the direction normal to the constant phase lines and cannot be near 
90^. The pressure fluctuation can be written more conveniently as 



( 100 ) 


where p' is the wide-band rms pressure fluctuation, and 7^ ^ dimensionless 

two-dimensional spectrum function. A two-dimensional spectrum function is 
needed to account for the distribution of energy through different orien- 
tations for the same frequency. What is being said here is that a harmonic 
coiUf/onent (cuAAl^>BA) of the freestream disturbance excites a Tollmien- 
Schlichting wave of the same frequency and orientation with an energy density 
proportional to the pressure-fluctuation energy density of the freestream 
component under the assumption that the imposed pressure fluctuation at the 

wall is the same as in the freestream. The proportionality factor is the 

2 

dimensionless frequency cdv/U^ , so that this line of reasoning has produced 
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the result that high frequencies are the most effective in producing insta- 
bility waves. There is no cross spectral transfer of energy. 


D. Effect of Freestream Turbulence on Transition 


1. Application of Amplitude Method 

In order to proceed further, we now restrict the freestream disturbance 

to turbulence, and furthermore consider the turbulence to be isotropic. As 

the pressure fluctuation appearing in (99) has been assumed to be the same as 

(92) 

in the freestream, isotropic turbulence theory gives 



oe 



( 101 ) 


Consequent ly , 


A 


2 


o 





cos^* 


( 102 ) 


As the X- Reynolds number does not appear in (102), we may take A^ to be the 
amplitude at the neutral-stability point. 

For a band of waves of the same frequency and different orientations, 
the two-dimensional (3 * 0) component will be the most amplified for fre- 
quencies along the envelope curve of f n(A/A^) , i.e., the maximum in(A/A^) at 
a given Reynolds number. These frequencies are given in figure 2 for the 
Falkner-Skan profiles. The calculation of the lateral bandwidth i[0(vxAJj^)^] 

in addition to fn(A/A ) , F and AF, involves the computation of four or 

o max max 

five times as many eigenvalues as are otherwise needed. Not enough numerical 
results are available to make a general statement about the behavior of the 
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lateral bandwidth, but for the flat-plate boundary layer it decreases slightly 
with increasing R. For example, when defined, as la AF, as the A0 at which 
has decreased to 1/e of its maximum value, it is 0.109 at R « 900 0.081 at 

R » 1600. Since this chaise has little influence on the transition Reynolds 
number, we will proceed on the basis that the two-dimensional component is 
the most amplified and that the lateral bandwidth is constant. Consequently, 
the substitution of (102) into (98) yields, with this simplification. 



(103) 


The three quantities on the second line ol (103) are all determined from 

9 

stability theory. The first two of these are needed for the e method, and 
the third is determined as explained in Section II-D.l without having to do 
any extra eigenvalue computations. 

The quantities in the first line of (103), except for R and the free 
constant C, are associated with the turbulence. It only remains to develop 
an expression for the spectrum function from isotropic turbulence theory. 

A convenient starting point is von Karman's normalized (to 2tt) interpolation 


Ej(kj^) = 4 



(104) 


for the one -dimensional spectrum function of the longitudinal velocity com- 
ponent. The dimensionless frequency twA/Uj^ has been replaced by the dimensionless 
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longitudltial wave number for the subsequent derivation. This procedure 
is permissible because all wave numbers have the convection velocity 
When Batchelor's theory ^ is applied to (104), a *7/3 rolloff is obtained 
for the one-dimensional pressure spectrum instead of -5/3 as in (104) . A 
normalized interpolation formula equivalent to von Karman's which has this 
behavior is 


Fj<kj) 



The one -dimensional spectrum is given in terms of the three-dimensional 
spectrum by 


k. 


k F^(k) dk , 


2 2 2 k 

where k = (kj^ ^ ^3 two-dimensional spectrum is given in 

terms of the three-dimensional spectnim by 


F2(k 


■/ 


P3(k) dk3 , 


2 2 k 

where kj^^ “ ^ ^ * 


F2(ki2 




dF, 


k (k^ - k 

*^12 ^ 12 ^ 


dk 


An approximate evaluation of (108), which uses (103) and is correct for 
kj^2 ^ ^ ^12 “ ® 


^2^^12^ = 1.78 [1 + (0.82 ' 


(105) 


(106) 


(107) 


(108) 


(109) 
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It may be observed that the rolloff in Is -10/3 as opposed to -7/3 for 
Fj^ and -5/3 for For use In (103) where ■ 8) is zero, kj ^2 replaced 

by oxA/Uj^. The two spectrum functions Fj^ and F 2 are shown in figure 8. It 
may be observed that the integral of F 2 is equal to Fj(0)/2 as required by 
the above equations . 

2 . Numerical Results 

All of the apparatus has now been assembled to make use of (103), but 
first the constant C must be evaluated. According to figure 7, Re^ = 2.8 x 10^ 
at T (or “ 0.001 for the flat-plate boundary layer, and with C = 1208 

A^ = 0.04 at this Reynolds number for = Uj^A/v “ 4 x 10 . In all subsequent 
calculations, C retains this value and transition is predicted to start whenever 
A^ first reaches 0.04. Any other value of A^ could have been chosen with a 
proportionate change in C. 

It is easy to see from (103) that for given values of R and R^, A^ 

2 2 

varies as T , and A^/T thus defines a single cur .. '^here are two effects 

of the scale Reynolds number R^ which tend to oppose each other. The first 

effect increases A^ through the appearance ''f R^ as the proportionality factor 
% 

relating the bandwidths referred to the boundary-layer scale length to the 

bandwidths referred to the turbulence length scale. Another way to look at 

this effect is that es the turbi'lence scale increases, the response curves 

spread over a greater req»;artcy range of the turbulence power spectrum, with 

the result that mor-s ereigy is included in the amplified band of frequencies 

and the transition Reynolds number is reduced. The seco? . effect is that as 

2 

R, increases, the unstable frequency band, which is fixed in terms of wv/U, 

(cf. figure 3), moves to Mgher dimensionless frequencies coAAJ^ ( * 

with a smaller energy density and the transition Reynolds number is increased. 


-56 



77-15 


/Qg\ 

The only reported measurements on the effect of scale on transition^ are 

for turbulence levels between 1% and 3%, and show a marked increase in the 

transition Reynolds nund>er with an increase in scale. The amplitude method 

4 

gives a decrease of 16% in Re^ at T » 0.02% as increases from 1 x 10 to 

4 

8 X 10 and the first effect prevails. At T * 1%» there is still a substantial 

linear amplification region before transition, but this region is not necessarily 

what is controlling transition. The Re^ of the amplitude method initially 

4 4 

decreases as R^ goes from 1 X 10 to 2 x 10 , but then increases by 16% as R^ 

4 

Increases further to 8 x 10 and the second effect prevails. This increasi 
is smaller, and Re^ higher, than in the experiment. 

Disturbance growth curves for the flat opiate boundary layer are shown 
in figure 9 for several turbulence levels. It is the steep slopes of these 
curves which are the basic Justification for the use of a method based on 
linear stability theory for transition prediction. Figure 10 shows Re^ as 
a function of T for six Falkner-Skan profiles, and the experimental flat-plate 
data are repeated from figure 7. The theoretical curve given by the amplitude 
method for 8*0 agrees well with these data for 0.09% <T < 0,27%. For 
T < 0.09%, transition is controlled by sound, and agreement with the theory 
is not expected. There is little point in being concerned about lack of agree- 
ment elsewhere, because the experimental data to properly test the theory do 
not exist at the present time. Only when enough data are available to make it 
possible to sort out the effects of T and R^ will it be possible to say whether 
the method given above properly accounts for the effect of freestream turbulence. 

The curves in figure 10 for the other values of 0 give some idea of the 
influence of a pressure gradient on transition. However, they cannot be used 
for the prediction of transition in a real boundary layer on the basis of local 
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values of 3» because Re^ depends on the amplification history of the disturbance 

Only if the particular value of 3 exists over a substantial portion of the 

boundary layer » a situation most likely to arise small favorable pressure 

gradients, can figure 10 be expected to give a quantitative estimate of the 

pressure gradient effect. Figure 11 is a crossplot of the results to show the 

effects of 3 for several turbulence levels. The effect of T on Re^ is almost 

independent of 3 for adverse pressure gradients, but for favorable pressure 

gradients T has a progressively reduced effect as 3 increases. 

It is of interest to know the value at transition of Jtn(k/k ) , i.e. 

o'max 

9 

the e factor n, according to the amplitude method. In the e method, n is of 

course equal to nine for all boundary layers and all turbulence levels; in the 
9 

modified e method n is given by (95) as a function of T for all boundary layers 

Only with the amplitude method is n a function of both T and the boundary layer. 

As 3 decreases from 3 ® 0 to -0.1 with T = 0.1%, n decreases from 8.2 to 7.7; 

as 3 increases to 0.2, n increases to 10.1. The effect of this variation in 

n on Re^ is more clearly seen if we refer to the dashed curve in figure 11 of 
9 

the e method. This curve is close to the T - 0.05% curve for 8 < -0.05, but 
with higher values of 3 it departs more and more until at 6 = 0.2 it corre^^ponds 
to T ~ 0.2%. Thus the amplitude method gives the result that turbulence levels 
which vary by a factor of four can all give n = 9 at transition depending on 
the value of 3* Only experiments on pressure gradient boundary layers with 
different turbulence levels can determine if the transition Reynolds number 
really varies with T and S as given by the amplitude method. 

The three jethods of transition prediction which have been presented in 
this section all use the same amount of computer time, because the eigenvalue 
computations are what determine the time. Only if the lateral bandwidth for a 
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two-dimensional boundary layer Is determined it each Reynolds number does 
the amplitude method require substantially more time than the other two 
methods. The amplitude method does have the disadvantage that a good 
deal must be known about the disturbance source, but this Is only a reflection 
of the physical situation that transition Is dependent on the type, Intensity 
and spectrum of the disturbance source, and not just on the boundary layer. 

The method has been set up so that as Information becomes available on the 
mechanisms by which Tollmlen-Schllchtlng waves are produced by external dis- 
turbances, more realistic luceractlon relations can be easily Incorporated 
into the computer program. It Is hoped that the development of a method which 
may for the first time offer the possibility of a rational prediction of trans- 
action will encourage the experimental work which Is so necessary to arrive at 
this long sought goal. 
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Table I. Operation of Eigenvalue Search Procedure 
e - 0, R = 1000, F = O.SxKT*. NSTEP - 80, = 8, A = 0.001 
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Table II. Properties of Falkner-Skan Boundary Layers 
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0.1372 

16,78 

0.473 

0,601 

0.579 

0.41 

0.2161 

17.65 

0.527 

0.2103 

-11.17 

0.513 

0.618 

0.633 


Table HI. Comparison of Spatial and Temporal Theories for Three Falkner-Skan 
Boundary Layers 
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♦ 

R 

F X 10* 


X 10® 

X 10® 

0® 

1600 

0.30 

0® 

-3.82 

-3.82 

45® 

2263 

0.15 

0® 

-3.17 

-3.17 




10® 

-3.12 

-3.07 




30® 

-3.32 

-2.88 




45® 

-3.82 

-2.70 

60® 

3200 

0.075 

0® 

-2 .46 

-2.46 




10® 

-2.42 

-2.39 




30® 

-2.59 

-2.24 




45® 

-2.98 

-2.11 




60® 

-3.82 

-1.91 

75® 

6182 

0.0201 

0® 

-1.37 

-1.37 




10® 

-1.36 

-1.34 




30® 

-1.49 

-1.29 




45® 

-1.75 

-1.24 




60® 

-2.32 

-1.16 




75® 

-3.82 

-0.99 


Table IV. Spatial Amplification Rates of Oblique Waves In the 
Flat-Plate Boundary Layer for Different Values of T. 
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77-15 



ure 3. Frequency-response curves of amplitude ratio 
at several Reynolds numbers for flat-plate 
boundary layer 



Figure 4. Reynolds number dependence of bandwidth of 
frequency-response curves for Falkner-Skan 
boundary layers 
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Figure 5. Direction of group velocity for oblique waves 
of constant wave number in flat-plate 
boundary layer, » 1600 



Figure 6* Comparison of theory with Schubauer-Skramstad 
measurements of the growth of six 
constant-frequency waves in flat-plate 
boundary layer 
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Figure 7. Effect of freestream turbulence on the transition Reynolds num- 
ber of the flat-plate boundary layer: o Schubauer-Skramstad ; 

□ Hall and Hislop^®^); A Wright and Bailey^®^); ODryden^^^); 
modified e^ method 
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Figure 8. One- and two-dimensional interpolation pressuve >ectra of isotropic 
turbulence. Fj^ la normalized to 2ir 
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Figure 9. Distcrbfln''e amplitude growth in flat-plate boundary layer 
according to amplitude meth'^d for several freestream tur- 
bulence levels, UjA/v « 4 - lO^ 
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'igure 10. Effect of freest! earn turbulence o-,i the transition 
Reynolds number of Falkner-Skan Loundary layers: 
UjA/v » i X 10^; o, a. A, O, experimental data 

(see caption of Figure 7 for sources); 

amplitude method 
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Figure 11. Effect of pressure-gradient parameter on the transition R€v- 
nolds number of Falkner-Skan boundary layers according to 
amplitude method for four levels of freestream tuibulcnce: 
UjA/v • 4 X 10^; amplitude method: nietnod 
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